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Abstract 

In this work a simple solid-on-solid (SOS) model of inhomogeneous films epitaxial 
growth is presented. The results of the computer simulations based on the random 
deposition (RD) enriched by a local relaxation which tends to maximize the num- 
ber of a particle-particle lateral bonds (PPLB) against barriers blocking diffusion 
are presented. The influence of strength of the atoms interactions and barriers for 
diffusion on film roughness and possible bridges between the layers in tri-layered 
sandwich system are considered. For magnetic layers with a nonmagnetic spacer the 
bridging is responsible for direct magnetic coupling between layers. Also variations 
of the spacer thickness, when scanning the film surface, is essential for uniformity of 
the magnetic coupling between the layers. Number of bridges for very thin spacer, 
and variation of the spacer thickness are discussed in this work. 
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1 Introduction 



In last decade the physics of growth thin solid films and ion/atoms behaviors 
on flat film surfaces became well understood and investigated due to rapid 
progress in experimental methods of preparing thin films (mainly MBE in 
UHV systems) and microscopy (mainly STM). Also many papers devoted to 
theory of growth of rough surfaces and interfaces and crystal were published 
(see [1-4] for review). 
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Here, we present a simple SOS model [5-7] of epitaxial growth adopted for 
inhomogeneous a/b/a-like systems. The film growth described as a two step 
mechanism. The first step is the randomly picked position of initial contact of 
the incident particle on the film surface. This is followed by second step which 
is a relaxation process of local migration leading to final position of the particle 
sticking to the film. We simplify the model to the limited surface diffusion 
and neglect migration to more distant sites then nearest neighbors (NN) only. 
After the relaxation process particle sticks for the rest of the simulation. Each 
particle is represented by a unit-volume cube which can occupy only discrete 
position in the lattice. The simple cubic symmetry is assumed. For a L x L 
substrate we deposit 9 a L 2 particles of kind a, followed by 9 b L 2 Mike particles 
and again 8 a L 2 of a particles. The nominal thicknesses of a and b layers, are 9 a 
and 9 b respectively. In relaxation process, each particle tends to maximize the 
number of PPLB. This tendency is slowed down by the barrier V for diffusion 
which decreases probability of atom movement. The above growth model may 
be implemented as the following flowchart: 

• for each newly arriving particle at random site r = 0, calculate the number 
of atomic pairs n r aa1 rf ab and rf bb at the place of the initial particle contact 
to the surface with all their four NN labeled by r = 1 ... 4, 

• calculate the particle total energies in all five positions (r = . . . 4): E r = 
rf aa E aa + rf ab E ab + n r bb E bb) where Eij is the bonding energy between i- and 
j-kind atoms, 

• evaluate probabilities p r oc exp(— E r ), r = ... 4, of picking out each of the 
five virtual final positions of the atom, 

• reduce the probability p r of movement into r = 1 . . . 4 by factor exp(V^), 
where V a and V b are the diffusion barriers for a- and Mdnd of atoms, re- 
spectively, 

• pick out one of five proposed sites for the atom with probability given by 
p r , (r = 0...4). 

Values of E and V are expressed in ksT units, where ks is the Boltzmann 
constant, and T denotes the absolute temperature. Diffusion barrier V is pos- 
itive, while the negative E are compatible with the assumed tendency of the 
system to maximize the number of PPLB. 



2 Results of simulation 

The simulations were performed on 500 x 500 large square lattice with periodic 
boundary conditions. The nominal thickness 9 a of both a-layers was set up 
to ten monolayers (ML), while average thickness of ft-spacer 6 b was varied 
from one to ten ML. The cases E aa = E ab = E bb — 0, or V a — Vb — > +oo 
which blocks any diffusion, yields the same results as RD model, with Poisson 
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distribution of surface/interface heights. 

2.1 Direct magnetic coupling 

Energy of magnetic coupling between magnetic layers separated by nonmag- 

— * — * 

netic spacer is often expressed by E = —KM 1 o M 2 energy term, where 

— * — * 

Mi and M 2 are magnetizations in magnetic layers. The coupling K depends 
strongly on nonmagnetic spacer thickness 9b and may follow exponential law 
K oc exp(— 9b/9o), or it may have an oscillatory character against 9b as it also 
is often observed. The interaction energy K(9b) manifests itself experimen- 
tally by modifications of some of magnetic properties such as susceptibility or 
ferromagnetic resonance. The coupling between magnetic layers separated by 
a nonmagnetic spacer was shown for example for NiFe/Cu/NiFe [8] and for 
Ni/Ag/NiFe samples [9]. In latter case for some samples also power decrease 
of K with 9b was observed. 

For simple RD model distribution of spacer thickness hb follows Poisson distri- 
bution P(h; 9) = 9 h /h\ ■ exp(— and probability of direct coupling between 
a-layers, corresponding to zero spacer height, decrease exponentially with av- 
erage spacer thickness 9 b : P(0; = exp(— 9b). However, in film preparation 
technology, the growth conditions seldom correspond to RD model and so we 
expect deviations from Poisson distribution as presented in Fig. 1. For differ- 
ent sets of model parameters for non-RD, we found that decrease of number 
of bridges (proportional to the direct coupling K) follow either exponential or 
the power law (Fig. 2). 

2.2 Spacer roughness 

Variation of the spacer thickness is also essential for uniformity of the mag- 
netic coupling between layers. We use root-mean-square a of surface heights 
as a measure of surface roughness. The dependence of aMayer roughness 
on roughness of each of its components is described by: a 2 ab = a 2 a + of + 
2{h a hb) — 2{h a )(hb)- For RD, successive film heights are uncorrelated and thus: 
a ab = °a + a b an d °~1 — 9{ for i = a, b, ab. The situation for non-RD case 
is more complex (see Fig. 3b-f). We would like to consider roughness ab de- 
pendence on the spacer thickness 9b for different model control parameters. 
For homogeneous films growth models roughness o scaling with film thick- 
ness with Family- Vicsek law [10]: o oc L a f(9/L z ) with f(x — > 0) oc x 13 and 
f(x — > oo) oc 1. For large enough lattice size, and not too large average film 
thickness, this dependence is given by a power law: a oc 9 13 . The dependence 
of exponent f3 on model control parameters was discussed in details in Ref. [6]. 
This way we are able to predict roughness cr a {9 a ) of the first a layer. However, 
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Fig. 1. Histogram of spacer heights for different model control parameters. Nominal 
thickness was 6 b = 4 ML. (a) E aa = -10, E ab = -10, E bb = -10, (b) E aa = 0, 
E ab = -10, E bb = 0, (c) E aa = 0, E ab = 0, E bb = -1, (d) E aa = 0, E ab = 0, 
E bb = -10, (e) E aa = -1, E ab = -10, E bb = -1, (f) E aa = 0, E ab = 0, E bb = 0. 

determination of a b is more difficult. Firstly, 6 b is small to guarantee clear 
dependence u b on 9 b . Secondly, the growth of b spacer is also controlled by 
E ab until Mike coverage is sufficient to creating 66-like PPLB when becomes 
dependent on E bb . Fig. 3 shows that a b increases with increasing of 9 b inde- 
pendently on E and V parameters. Deviations from basic power scaling law, 
however, may be observed. 



3 Conclusion 



We found from computer simulations that distribution of spacer heights around 
average thickness changes from Poisson distribution (or bell-shaped for larger 
thickness), to peak-shaped with a decrease of diffusion barriers V and/or in- 
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Fig. 2. Exponential (a) and power-law-like (b) dependence of the probability of 
bridging on the spacer thickness 9b for different model control parameters E aa , E a b 
and shown in rigth upper corners. 

creasing E, for a tendency of particles to create more PPLB. The latter case 
helps to produce more uniform of magnetic coupling between layers in the 
sandwich tri-layer systems. 

Energy of magnetic coupling between two magnetic moments fa and fa may 
be expressed by Heisenberg term E = —Kfa o fa, where K is the coupling 
constant. It is common to assume that this interaction is a short range one. 
For direct exchange interaction J, K oc J and J ^ only if hb = 0. The direct 
coupling means that hf, = and strength of coupling K between magnetic 

— * _ — * 

layers with magnetizations Mi = J2fa an d M 2 = J2fa ma Y by evaluated as 
K = NJ, where iV = P(0; 0&) • L 2 is the number of bridges between magnetic 
layers. Usually, the decrease of K with increasing spacer thickness follows ex- 
ponential or power law. We found from computer simulations that number of 
bridges responsible for direct coupling is compatible with the above predic- 
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Fig. 3. The lines of +, x, * and □ correspond to a a , a b , a ab , and a aba , respectively, 
(a) E aa = 0, E ab = 0, E bb = — a RD model, (b) E aa = 0, E ab = -1, £ b6 = -10, 
(c) £ aa = -1, E ab = 0, £ 66 = -1, (d) E aa = -1, ^ab = -10, £ 66 = 0, (e) E aa = -1, 
E ab = -10, £ b6 = -10, (f) E aa = -10, E ab = 0, £ b6 = 0. 

tions, and either exponential or power law may be obtained for specific set of 
model parameters. 



Acknowledgements 



This work and machine time in ACC-CYFRONET-AGH is financed by Polish 
Committee for Scientific Research (KBN) with grants no. 8 T11F 02616 and 
KBN/S2000/AGH/069/1998, respectively. 



References 



6 



[1] J.-F. Gouyet, M. Rosso, B. Sapoval, in Fractals and Disordered Systems, A. 
Bunde, S. Havlin, eds. (Springer- Verlag, Berlin, Heidelberg, 1991) p. 229. 

[2] H. J. Herrmann, Phys. Rep. 136 (1986) 153. 

[3] P. Meakin, Phys. Rep. 235 (1993) 189. 

[4] A. C. Levi, M. Kotrla, J. Phys.: Condens. Matter 9 (1997) 299. 

[5] A. Z. Maksymowicz, M. S. Magdori, J. S. S. Whiting, Comp. Phys. Comm. 97 
(1996) 101. 

[6] K. Malarz, A. Z. Maksymowicz — |cond-mat/9905110| {Int. J. Mod. Phys. CIO 

- in print). 

[7] K. Malarz, A. Z. Maksymowicz — |cond-mat/9905111| (Int. J. Mod. Phys. CIO 

- in print). 

[8] A. Z. Maksymowicz, J. S. S. Whiting, M. L. Watson, A. Chambers, J. Magn. 
Magn. Mat. 94 (1991) 109. 

[9] A. Layadi, J. O. Artman, R. A. Hoffman, C. L. Jensen, D. A. Saunders, B. O. 
Hall, J. Appl. Phys. 67 (1990) 4451. 

[10] F. Family, T. Vicsek, J. Phys. A18 (1985) L75. 



7 



